Strong polaritonic interaction between flux-flow and phonon resonances in Bi2Sr2CaCu208+x 
intrinsic Josephson junctions: Angular dependence and the alignment procedure. 
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Bi2Sr2CaCu208+jc single crystals represent natural stacks of atomic scale intrinsic Josephson junctions, 
formed between metallic Cu02-Ca-Cu02 and ionic insulating SrO-2BiO-SrO layers. Electrostriction effect 
in the insulating layers leads to excitation of c-axis phonons by the ac-Josephson effect. Here we study exper- 
imentally the interplay between and velocity matching (Eck) electromagnetic resonances in the flux-flow state 
of small mesa structures with c-axis optical phonons. A very strong interaction is reported, which leads to 
formation of phonon-polaritons with infrared and Raman-active transverse optical phonons. A special focus in 
this work is made on analysis of the angular dependence of the resonances. We describe an accurate sample 
alignment procedure that prevents intrusion of Abrikosov vortices in fields up to 17 Tesla, which is essential for 
achieving high-quality resonances at record high frequencies up to 13 THz. 

PACS numbers: 74.72.Hs, 74.78.Fk, 74.50.+r, 85.25.Cp 



Intrinsic Josephson junctions (IJJs) |1 1 in high temperature 
superconductors are considered as possible candidates for re- 
alization of high power THz-frequency radiation sources (2]- 
[7 |. IJJs can provide several advantages: (i) A large number of 
IJJs can easily be integrated. A strong superradiant emission 
can take place if all IJJs are phase-locked in the symmetric 
in-phase mode, (ii) Due to the large energy gap in Bi-2212 
A ~ 30 - 40meV |8|, IJJs can comfortably operate in the 
frequency range up to ~ 20 THz. (iii) Stacking of IJJs may 
lead to cascade amplification of non-equilibrium population, 
needed for achieving population inversion and lasing |4 |. (iv) 
A large variety of emission mechanisms are available. The 
ac-Josephson effect facilitates electromagnetic wave emission 
via fluxon-induced self-oscillations at geometrical resonances 
(61 m, or via the Josephson flux-flow emission ||7l[T0l. Be- 
sides, coherent phonon emission can take place either without 
the ac-Josephson effect via recombination of non-equilibrium 
quasiparticles |4|, or via the ac-Josephson effect due to the 
electrostriction phenomenon Llll . 

In this work we investigate the mutual interaction between 
phonons and electromagnetic waves in the flux-flow state of 
small Bi-2212 mesa structures. The ac-Josephson effect leads 
to excitation of c-axis optical phonons in Bi-2212 single crys- 
tal by means of electrostriction: Ions in the polar insulating 
SrO-2BiO-SrO layers, where the oscillating electric field is 
concentrated, are shaken at the Josephson frequency. Both Ra- 
man and infrared (IR) - active c-axis phonons can be excited 
lfTTUT3l . The symmetry of the phonon is directly related to 
the symmetry of the electromagnetic mode: Raman phonons 
are excited by the out-of-phase mode, while infrared- active 
phonons are excited by the in-phase mode. The strongest in- 
teraction between electromagnetic waves and phonons occurs 
at the transverse optical phonon frequency. This leads to for- 
mation of phonon-polaritons with very slow group velocities. 
We argue that resonances with IR-active phonons can stabilize 
the in-phase electromagnetic mode in IJJs, needed for achiev- 
ing coherent superradiant flux-flow emission from Bi-2212 
mesas. The appearance of phonon-polaritons also demon- 



strates the existence of an unscreened polar response, which 
may lead to the strong electron-phonon coupling in cuprates. 

We study small mesa structures, made at the surface of 
freshly cleaved Bi-2212 single crystals. Mesas are made 
by photolithography, argon milling and focused ion beam 
trimming down to sub-micron sizes. Two types of crys- 
tals are studied: slightly underdoped pure Bi-2212 crys- 
tals with Tc - 82 K and slightly overdoped lead-doped 
Bi2-3;PbySr2CaCu208+x (Bi(Pb)-2212) crystals with ^ 
90 K. The most noticeable diflTerence between those crys- 
tals is in the critical current density, which is an order of 
magnitude larger for Bi(Pb)-2212 (/c ~ 10^ A/cm^). Since 
the amplitudes of the phonon resonance, reported below, are 
proportional to JI |14|, they are much more pronounced in 
Bi(Pb)-2212 crystals |11|. In order to avoid repetitions we 
address the reader to our previous works, in which we stud- 
ied properties of similar mesas in intense magnetic fields, in- 
cluding: the Fraunhofer modulation of the Josephson current 
115], velocity matching (Eck), geometrical (Fiske) |10| and 
phonon-polariton |11| resonances in the flux-flow state, their 
temperature dependence |[T6ll and the intrinsic tunneling mag- 
netoresistance of IJJs |17|. Below we will focus on the angu- 
lar dependence and describe the procedure for accurate sam- 
ple alignment, which is a prerequisite for observation of high 
quality resonances. 

Figure [T] represents the angular dependence of the current- 
voltage characteristic of the Bi(Pb)-2212 mesa P22/m5b with 
the area ~ 1.3 x 0.9yum^, containing N ^ \\ IJJs. For 
the field parallel to the aZ^-plane, (a) a = 0, the quasipar- 
ticle (QP) branches are completely suppressed and a pro- 
nounced flux-flow step appears at V ^ 0.1 V. It corresponds 
to the Eck resonance, when the speed of propagating Joseph- 
son vortices (fluxons) is equal to the speed of the slowest 
out-of-phase electromagnetic wave |[T0l[TT][T6l. Strong hys- 
teretic resonances above the Eck step at higher voltages, 0.1 V 
< y < 0.2 V, are due to excitation of c-axis optical phonons 
|[TTHT4l [TSl . Already at a slight misalignment of the field by 
1-2 degrees (b, c) the resonances are strongly dampened due to 
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FIG. 1. (Color online). Angular dependence of current- voltage characteristic for a Bi(Pb)-2212 mesa P22/m5b at 7 = 1.7 K and H - 17 T. 



penetration of Abrikosov vortices. The amount of Abrikosov 
vortices increases with further increase of the angle (d), they 
effectively pin Josephson fluxons, which leads to a significant 
increase of the critical current, suppression of the flux-flow 
and a reappearance of QP branches (d, e). With further in- 
crease of the angle the amount of Abrikosov vortices becomes 
so large that their normal cores start to suppress the supercon- 
ducting order parameter. This results in the pronounced nega- 
tive magnetoresistance 1 17 1 and the suppression of the critical 
current (f-h). 

Figure [2] represents difl'erential conductance curves dl/dV 



0.004 




V/N (mV) 



FIG. 2. (Color online). Tunneling conductance versus voltage per 
junction for a Bi(Pb)-2212 mesa P22/m4a at T = 20 K and // = 16 T 
in the ab-plstnc and in the c-axis direction. 



versus voltage per junction V/N for another Bi(Pb)-2212 mesa 
(P22/m4a 1.2 x 2.0 jim^, N ^ 55) for field along the a/?-plane 
and in the c-axis direction. Phonon resonances are seen as 
peaks in dl/dV for a field along the a/?-plane. They com- 
pletely disappear for field in the c-axis direction, clearly show- 
ing the importance of sample alignment. 

Table |T| summarizes all observed phonon resonances for the 
Bi(Pb)-2212 crystal. In total fourteen phonon resonances with 
varying amplitudes could be distinguished. This is consis- 
tent with the expected number of c-axis optical phonons in 
the tetragonal Bi-2212 (seven IR and seven Raman) |[T9ll . In- 
deed, by comparison with optical spectroscopy [191, most of 
the observed phonons are easily identified with the known 
Raman and IR-active c-axis optical phonons lfTTHT3ll . Few 
unidentified phonons may be disorder induced, due to the or- 
thorhombic distortion, or the incommensurate superstructure 
in the crystal |[T9ll . 

From Fig. [T] it is obvious that an accurate sample align- 
ment is essential for the study of phonon resonances in intense 
fields. Therefore we want to describe the developed unam- 
biguous alignment procedure in more details. Figure [3] shows 
the angular dependence of the dc-resistances R^c = V/I, mea- 
sured at diff'erent probe currents / (a) low / = 0.05 mA, (b) 
moderate / = 0.1 mA and (c) high 1 = 3 mA for the Bi(Pb)- 
2212 mesa P22/m3 at T ^ 4 K and at three diff'erent magnetic 
fields. The angular dependence of the magnetoresistance can 
be understood from I-Vs in Fig. [T] The sharp maximum at 
a = is caused by the flux-flow contribution. The flux-flow 
voltage is rapidly suppressed at small a misalignment a - 1° 
due to the pinning by Abrikosov vortices. At larger angles the 
magnetoresistance strongly depends on the bias: at a low bias 
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assignment Bi:CuCaSr CuSr SrCu Sr:Cu "disorder" Ca:Sr - - "disorder" 03:01 - 01:Ca03 01:Sr 

^ very weak 

^ disorder-induced IR-active or Big Raman-active 

^ disorder-induced IR-active or Ai^ (03) Raman-active 

TABLE I. Experimental voltage positions of phonon resonances. The mean values for all mesas at a Bi(Pb)-2212 crystal P22 are given. Also 
given are types (Raman-active resp. IR-active), symmetries and assignments of phonons. 
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FIG. 3. (Color online). Angular dependence of dc-resistances, measured at different probe currents / (a) low / = 0.05 mA, (b) moderate 
/ = 0. 1 mA and (c) high / = 3 mA for the Bi(Pb)-2212 mesa P22/m3 at T 4 K and at three different magnetic fields 4 T (sohd), 14 T (dashed) 
and 17 T (dashed-dotted line). 



the resistance increases due to the suppression of the critical 
current, but at a higher bias it is monotonously decreasing due 
to the negative QP magnetoresistance lfT7ll2T1l . 

Even though the sharp flux-flow peak at low bias seems to 
be advantageous for the sample alignment, we found that such 
an alignment is unreliable because of the field lock-in along 
the a/?-plane, which may lead to a hysteresis up to a few de- 
grees |20|. In Fig.|3](a) it can be seen as a slight asymmetry 
of the minima in the /?dc(Q^) curve at // = 17 T. With in- 
creasing field the angular range of the hysteresis decreases. 
At H = 14T it is only ~ 0.2°. However this corresponds to 
almost 500 G c-axis field that is suflftcient for penetration of 
Abrikosov vortices. Therefore the low bias alignment is not 
good enough. On the other hand, we observed that the align- 
ment at high bias, as shown in Fig. [3] (c), is working much 
better because it is free from the hysteresis. This is caused 
by an eflfective depinnig of Abrikosov (pancake) vortices as 
a combination of self-heating that leads to thermal-assisted 
flux creep, a gradient force due to perpendicular magnetiza- 
tion in pancake vortices with respect to applied field and a 
drag force from Josephson vortices. Furthermore, the max- 
ima Rdcio:) become significantly sharper at high bias. For the 
curves ai H = 14 T the full width at half maximum is 3.4° 
at low bias. Fig. [3] (a), compared to only 1.5° at high bias, 
Fig.|3](c). Therefore we use high bias magnetoresistance for 
the accurate sample alignment. 



Figure |4] shows the measured critical current vs the in-plane 
magnetic field for (a,b) a moderate size Bi-2212 mesa and (c) 
a small Bi(Pb)-2212 mesa P22/m4a. A clear Fraunhofer mod- 
ulation is seen in Figs.|4](a) and (c). A transition from the half 
flux quantum to the flux quantum Oq periodicity occurs with 
increasing field due to the reconstruction of the Josephson 
fluxon lattice from the triangular to the rectangular 03112211 . 
Panel (b) shows IdH) for the same mesa and alignment as in 
(a) after subjecting the sample to a larger field ~ 8 T. Appar- 
ently, the Fraunhofer modulation is completely destroyed by 
the intrusion of Abrikosov vortices. In Fig. [4] (c) the Fraun- 
hofer modulation persists up to 13 T, demonstrating that with 
our rigorous alignment procedure we can prevent Abrikosov 
vortex intrusion even at very high fields. 

Figure [5] shows the I-V oscillogram for the Bi(Pb)-2212 
mesa P22/m5b at T = 1.6 K in the in-plane magnetic field 
of lOT. Here a flux-flow branch ending with a nearly verti- 
cal Eck resonance is clearly seen. Above it a large variety 
of phonon resonances appears. Those resonances are divided 
into groups ofN= 11 branches due to one-by-one switching 
of IJJs into the resonant state. For clarity we painted different 
families of branches by different colors. 

Figure |6] (a) represents the magnetic field dependencies of 
collective resonance voltages (all junctions at the same reso- 
nance) for one of the studied Bi(Pb)-2212 mesas. Positions of 
all fourteen phonon resonances are indicated by dotted lines. 
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FIG. 4. (Color online). Fraunhofer patterns /c vs. in-plane H for a moderate size Bi-2212 mesa with the length L - SAjdm dX T - 
1.6 mathrmK: (a) For field sweeps from zero to 2.1 T and backwards, (b) After trapping Abrikosov vortices at higher field due to slight 
misalignment. Diff'erent curves in (a,b) represent separate runs, (c) For a small Bi(Pb)-2212 mesa L ^ 1.2 jum after rigorous alignment at large 
bias. The clear Fraunhofer modulation indicates the lack of Abrikosov vortex intrusion up to 13 T. Arrows indicate the field for introduction 
of one flux quantum Oq in a single junction. A transition from Oo/2 modulation at low fields to Oq at large fields is due to the reconstruction 
of the Josephson fluxon lattice from triangular to rectangular. 
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FIG. 5. (Color online) I-V oscillogram for the Bi(Pb)-2212 mesa 
P22/m5b 3i H = 10 T aligned along the aZ^-plane. Diff'erent fami- 
lies of phonon-polariton resonances are painted in difi'erent colors. 
They are due to one-by-one joining of junctions to the resonant state. 
The extraordinary large amplitude and hysteresis of the resonances 
indicates their high quality factors and oscillation amplitudes. 

As discussed in Ref. flT|, the V(H) plot can be translated 
into the dispersion relation with the frequency given by the 
ac-Josephson relation and the wave length by the periodicity 
of the Josephson fluxon lattice: 

/ = V/A^Oo, (1) 
k = 2nHs/^o, (2) 

where s = 1.5 nm is the stacking periodicity of Bi-2212. 
Therefore the V-H plot represents the dispersion relation of 
electromagnetic waves in the mesa. The lowest branch (black) 
represents the Eck resonance. Initially it increases linearly 
with field, corresponding to the constant speed of electromag- 



netic waves ^ 4 x lO^m/s, close to the out-of-phase speed 
ifTOl . However, at // > 12T it saturates upon approaching the 
second phonon mode. At this point the group velocity of elec- 
tromagnetic waves becomes very slow. This is a well known 
signature for formation of phonon-polaritons with the trans- 
verse optical phonon 1 18|. 

The polaronic dispersion relation is determined by the di- 
electric function, which is influenced by the interaction of 
propagating electromagnetic waves in the dielectric barrier 
with optical phonons: 



j ^Toj - ^yj^ 

Here ojtoj are frequencies of transverse optical phonons, S j 
and Yj are the corresponding oscillator strength and the damp- 
ing parameter, respectively. 

Figure [6] (c) shows a simulated dielectric function with 
characteristic phonon frequencies aligned to the positions of 
the fourteen observed resonances. The oscillator strengths 
were chosen so that s(0) = + is equal to the exper- 

imental value ~ 10 and damping parameters were chosen so 
that the corresponding dispersion relation oj = c(0)k/ ^Js{cjS), 
shown in Fig.|6](b) is looking alike the experimental data in 
Fig. [6] (a). Figure |6](b) shows the obtained dispersion rela- 
tion for positive group velocities with scales in frequency and 
the wave number corresponding to voltage and field scales in 
Fig.|6](a) according to Eqs. ([T]) and (|2]). Even though Fig.[6](b) 
is not meant to be a direct fit to Fig.]6|(a), it demonstrates that 
peculiarities of experimentally observed V{H) characteristics 
can be fairly well explained by the simple expression Eq. ([3|. 
In particular, the saturation of V{H) and small group veloc- 
ities are caused by the divergence of the dielectric function 
when the ac-Josephson frequency is approaching any of the 
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FIG. 6. (Color online), (a) Magnetic field dependence of voltages (per junction) of diff'erent resonances for one of the Bi(Pb)-2212 mesas. 
It can be translated into the dispersion relation f(k) according to Eqs. ([TJ and ([2|, which correspond to the scales in panel (b). Positions of 
fourteen characteristic phonon frequencies are marked by dotted lines. Panels (b) and (c) represent the calculated dispersion relation and the 
corresponding polaronic dielectric function, given by Eq. ([3]). 



active transverse-optical phonon frequencies. Also an abrupt 
jump-like cancellation of resonances at higher voltages (see 
Fig. [5]) is caused by the negative value of s, which prevents 
propagation of electromagnetic waves with such frequencies 
and leads to switching of I-Vs from the resonant flux-flow to 
the non-resonant quasiparticle state. 

In conclusion, we have described the accurate alignment 
procedure and the angular dependence of phonon-polariton 
resonances between electromagnetic waves, generated in the 
flux-flow state, and c-axis optical phonons in Bi-2212 single 
crystals. The phenomenon is caused by electrostriction in the 
insulating SrO-2BiO-SrO layers. Our work has several im- 
plications. For realization of the coherent flux-flow oscillator 
it is necessary to stabilize the in-phase state with the rectan- 
gular Josephson fluxon lattice, which is usually unstable be- 
cause of fluxon-fluxon repulsion. Strong interaction with in- 
frared optical phonons can stabilize the in-phase state. The 
stabilizing role of phonons is seen already from the fact that 
we observed the flux-flow phenomenon up to the record high 
frequencies ~ 13THz, see Fig. |6](a) and (b). Formation of 
phonon-polaritons demonstrates also that a residual polar re- 
sponse is left in cuprates even in the superconducting state. 
This demonstrates that Bi-2212 single crystals represent nat- 
ural atomic metamaterials composed of two-dimensional su- 
perconducting Cu02-Ca-Cu02 layers stacked with polar insu- 
lating SrO-2BiO-SrO. This can lead to an extraordinary strong 
electron-phonon coupling due to the combination of a long 
range Coulomb interaction in 2D metallic planes with an un- 
screened polar response in the adjacent polar insulator lay- 
ers. This may have significance for understanding the pairing 
mechanism in cuprate superconductors. 
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